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Abstract. We discuss the (right) eigenvalue equation for H, C and K linear quaternionic 
operators. The possibility to introduce an isomorphism between these operators and 
real/complex matrices allows to translate the quaternionic problem into an equivalent 
real or complex counterpart. Interesting applications are found in solving differential 
equations within quaternionic formulations of quantum mechanics. 
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I. INTRODUCTION 

The full understanding of the subtleties of the quaternionic eigenvalue problem still represents an in¬ 
triguing challenge for mathematicians and physicists. The recent study of the eigenvalue problem for 
complex linear quaternionic operators Q played a fundamental role in solving quaternionic differential 
equations Q. In the last few years, interesting applications of quaternionic analysis and linear algebra 
were investigated in quantum mechanics ||^. In particular, the solution of the Schrodinger equation 
in presence of quaternionic perturbations was explicitly given for constant potentials and deviations 
from standard (complex) quantum mechanics discussed [^. In this paper, we aim to complete the 
study begun in ref. where preliminary steps in solving the eigenvalue problem for complex linear 
quaternionic operators were traced. In order to extend to the R-linear case the results obtained for the 
H and C-linear quaternionic matrices, we have to introduce a system of coupled equations which rep¬ 
resents the new eigenvalue problem for R linear quaternionic operators. It is important to observe that 
no attempt to develop a complete theory of the quaternionic eigenvalue problem has been made here, 
this exceeds the scope of our paper. A satisfactory discussion of the eigenvalue problem for quater¬ 
nionic operators is at present far from being given. We could have directly investigated the eigenvalue 
equation in the quaternionic space, but we have preferred a more practical approach and chosen to 
handle the problem by finding a more familiar real or complex space isomorphic to the quaternionic 
one. We shall show that the isomorphism between H, C, and M-linear quaternionic operators and 
real/complex matrices immediately allows to translate the quaternionic (right) eigenvalue problem 
in a corresponding real or complex counterpart. The study of the new translated problem gives im¬ 
portant information about the quaternionic solution. The results obtained are very useful in solving 
polynomial and differential equations. This could represents a fundamental step in understanding the 
potentiality of using quaternions in formulating quantum mechanics (by investigating quaternionic 
deviations from the standard theory ||^,|| ) and gauge theory (by suggesting new unification groups § ). 

Throughout the paper we shall denote by R, C and H the sets of real, complex and quaternionic 
numbers, R C C C H, and by U[n,X] and M[n, X], respectively, the n-tuples and the nxn matrices 
over X. Linear quaternionic operators will be distinguished by their linearity from the right. In what 
follows, the notation Ox stands for quaternionic operators linear (from the right) over the field X. 
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II. QUATERNIONIC ALGEBRA AND LINEAR OPERATORS 


We now introduce the quaternionic algebra and some useful properties of H, C and R-linear operators. 
The (real) quaternionic skew-field HI is an associative (division) algebra of rank 4 over K, 

q = qo + iqi+jq2 + kq:i , 90 , 1 . 2 , 36 ®, ( 1 ) 

where 

*2 = f = e = ^Jk = -1 , ( 2 ) 

endowed with an involutory anti-automorphism {conjugation) 


q 9 = 9o - * 9i - j 92 - fc 93 ■ 


Due to the non-commutative nature of quaternions, we must distinguish between the left and right 
action of the quaternionic imaginary units i, j and k. To do it, we introduce the operators 

k'fx — (-1, , Lj^ Lj^) and Rfj, — (1, Ri , Rj , R ]^) , /i = 0,l,2,3, 

which act on quaternionic vectors ip G y[n,]HI] in the following way 

L^ip = h^ip and R^ip = iph^, = {l,i, j,k) . 


These operators satisfy 


Li = L] = Ll= L,L,Lk = Rl = R] = Rl = RkRjR^ = -1 , (3) 

and 

[L^,i?,]=0, /i,i/ = 0,l,2,3 . (4) 

Note that H-linear quaternionic operators acting on a finite n-dimensional quaternionic vector space, 

Oa {ipiqi + '0292) = (Oh'0i) 91 + (Oh02) 92 , 91,2 G HI , 01,2 6 V[n, H] , 

are in one to one correspondence with n x n quaternionic matrices: 

Oa Ma 6 ikf[n,IHI] . (5) 

Consequently, ffi. and C-linear quaternionic operators ^ 


Oa {tpiri + 1 p 2 r 2 ) = (C’r0i) D + (C’e02) ?’2 , D,2 G K , 01,2 G V[n, H] , 

Oc (0lCl -p 02C2) = (C’c0l) Cl -p (Oc02) C2 , Cl,2 G C , 01,2 G V[n, H] , 

can be represented by n x n quaternionic matrices Ma and right acting operators as follows 

3 1 

Oa ^ Ma = '^ M^,h R^ and Oc ^ Me = ^ Ms,„ Rs . (6) 

/J.—0 s—0 


Thus, R-linear quaternionic operators consist of right multiplication by quaternionic numbers (/i = 
0,1,2, 3) whereas C-linear quaternionic operators are restricted to right multiplication by complex 
numbers (s = 0, l). 


III. THE EIGENVALUE PROBLEM 

In this section, we briefly discuss the left and right eigenvalue equation for H, C and K.-linear quater¬ 
nionic operators. As explicitly shown below, the conceptual difficulties which characterize the left 
eigenvalue problem readily disappear by resorting to right eigenvalues. The need to apply similarity 
transformations on H, C and M linear quaternionic matrices introduces complex or real constraints 
on the right eigenvalues. The choice of complex or real (right) eigenvalues will be extremely useful in 
finding a practical method of resolution and manipulating quaternionic matrices. 
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lll-A. LEFT EIGENVALUES 

The left eigenvalue problem for H linear quaternionic operators reads 

Oui’ = qip, Ip € V[n,M\ , q € U . (7) 

This problem has been recently studied in the mathematical literature Nevertheless, no sys¬ 

tematic way to approach the problem has been given. We point out some difficulties which appear in 
solving the left eigenvalue equation. 

• Similarity transformations. In finding the solution of Eq. (I), a first difficulty is represented by the 

impossibility to apply similarity transformations, G without losing the formal structure 

of the left eigenvalue equation. In fact, by observing that Saq ^ q Sji, the quaternionic matrices 

Mb and (Sh M^ (8) 

do not necessarily satisfy the same eigenvalue equation. Consequently, we can have quaternionic ma¬ 
trices with the same left eigenvalue spectrum, but no similarity transformation relating them. Explicit 
example are found in ref. M. 

• Hermitian operators. Let Oh be an hermitian quaternionic operator, ip the eigenvector corresponding 
to the eigenvalue q. By using Eq. 0 and denoting by (iflip) the inner product in E[n,]HI], we obtain 

0 = {OB-tp\-tp) - i'iplOBtp) => 0 ^ {qip\tp) - {tplqtp) {q-q) {tp\tp) . 

Consequently, the left eigenvalue problem for hermitian operators could admit quaternionic solu¬ 
tions p. 

• Square operators and eigenvalues. As a last difficulty in the use of left eigenvalues, we observe that if 
Ip is an Ob eigenvector with eigenvalue g, it will not necessarily be an O^ eigenvector with eigenvalue 
q^. In fact, 

Olip = O^qip q Gulp = q^'ip ■ 

lll-B. RIGHT EIGENVALUES 

The right eigenvalue equation for H-linear transformations reads 

C>h'0 = '<pq , V' S H] , g G H . (9) 

Such an equation can be reduced to a right complex eigenvalue equation rephasing the quaternionic 
eigenvalues by unitary quaternions u, 

Ob^Pu = ipuuqu = pjaz , z G C . (10) 

This trick obviously fails for complex and real linear transformations. In fact, due to the presence of 
the operators Ri in Oc and R = {Ri, Rj, Rk) in Or, we cannot apply unitary transformations from 
the right. Observe that 

(0,,cV') u ^ Ob,c{iPu) , u G H . (11) 

The failure of the associativity in Eq. & suggests that we should consider complex eigenvalue equa¬ 
tions for C-linear quaternionic operators, 

Oc'ip = p’z , zGC, (12) 

and real eigenvalues for R-linear quaternionic operators, 

Oxtp = tpr , r G K. . (13) 

These equations are formally invariant under C and R-linear similarity transformations. Moreover, it 
can easily be proved that 

It is important to note here that R-linear quaternionic operators admit real eigenvalues only in par¬ 
ticular cases. Thus, Eq.(|^ has to be generalized. As shown later, a satisfactory discussion of the 
eigenvalue problem for K-linear quaternionic operators will require the use of a system of coupled 
equations. 
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IV. CANONICAL FORMS 

In this section, following the procedure introduced in the paper of ref. Q, we discuss the canonical 
forms for H, C and R-linear quaternionic matrices. The results we will establish find an immediate 
application in the theory of quaternionic differential operators. In fact, by using the canonical form 
Jx of a given matrix we can readily obtain the exponential 

exp [Mx x] = Sx exp [J* x] S~^ 

and consequently, avoiding tedious calculations, to solve quaternionic differential equations with con¬ 
stant coefficients |^. 

IV-A. H-LINEAR MATRICES 

While matrices over commutative rings have gained much attention, the literature on matrices with 
quaternionic entries is often fragmentary. The main difficulty is that, due to the noncommutative na¬ 
ture of quaternions, the standard method of resolution breaks down. Consequently, finding eigenvalues 
and canonical forms represents a more delicate problem. The recent renewed interest in quaternionic 
matrix theory 00 and its applications shed new light on this intriguing research field. To 
facilitate access to the individual topics, we recall the main properties of H-linear quaternionic ma¬ 
trices and repeat the relevant theorems from without proofs, thus making 

our exposition self-contained. 

In approaching the problem of diagonalization we have to consider a right quaternionic eigenvalue 
equation. In fact, from 


Mailjk = ipkQk , fc = l,...,n, tpk &V[n,M] , (14) 

in the case Mj, is diagonalizable, we immediately get the following matrix equation 

V'2, = Sa[tpl,tp2,—,i’n]DB , 

where 

Hh = diag[gi,92,---,<7n] 

and ..., V'n] is defined by ColkiSn) = ipk- 

Consequently, the diagonalization of the matrix Mh 

= SjnDtiSfj 

is obtained by solving the corresponding right eigenvalue problem. It is important to note here that 
we have infinite ways to diagonalize a quaternionic matrix Mj,, 

diag[ui ,U 2 ,...,Un] diag[( 7 i, 92 , • ■ •, <7n] diag[ui, 122 , • ■•,««] • 

Geometrically speaking this means that Im[g] can arbitrarily be fixed on the sphere of ray |lm[( 7 ]|. By 
a particular choice of the unitary matrix 

Ua = diag[ui,U 2 , ...,Un] 

we can set a preferred space direction, for example the positive i axis, and consequently a complex 
(positive) eigenvalue spectrum. 

Let us now briefly recall some properties of the eigenvalue spectrum of H linear quaternionic 
matrices. By using the symplectic decomposition of the matrix Mj, 

Ma = Ml + jM 2 , Ml 2 £ M[n,C] , 

and the symplectic decomposition of the vector if, 


'tf = 'tfi+ jif 2 , '01,2 G V[n, C] , 
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we can rewrite Eq.(^) in the following (complex) form 

Mb'!/' = , (15) 

where 

(m 2 ~m\) ^M[2n,C] and (^6) 

The mapping 

(17) 

is an isomorphism of the ring of quaternionic matrices Mb into the ring of the corresponding com¬ 
plex counterparts Mb. It is important to observe that this isomorphism do not preserve the inner 
product of eigenvectors |^. Nevertheless, the choice of a complex projection of quaternionic inner 
products |l^] opens the door to interesting applications in relativistic quantum mechanics ||^. The 
complex orthogonality of quaternionic eigenvectors (instead of a quaternionic orthogonality) implies 
a doubling of solution in the two-dimensional quaternionic Dirac equation p^ . The four (complex) 
orthogonal quaternionic solutions describe particle/antiparticle with spin up/down. The use of com¬ 
plex inner products is also a fundamental ingredient in the formulation of gauge theories by geometric 
algebras [||. 

The next theorem states the main property of the eigenvalue spectrum of H linear quaternionic 
matrices, for a detailed discussion see refs. |p4|,p^,p^. 

Theorem 1 . 

Let Mjh be the matrix given in Eq. Then, its eigenvalues appear in complex conjugate pairs. 

By using the result of Theorem 1 and the Gram-Schmidt method, we can readily obtain the triangular 
form for H-linear quaternionic matrices. 

Theorem 2 . 

Every M,, is unitarily similar to an upper triangular matrix. 


Moreover, a Jordan form can be given for every H-linear quaternionic matrix. 


Theorem 3 . 

Every n x n matrix with real quaternion elements is similar under a matrix transformation with real 
quaternion elements to a matrix in (complex) Jordan normal form with diagonal elements in the com¬ 
plex field. 


To prove Theorem 3, we can use the isomorphism defined in & To any corresponds a 2n x 2n 
matrix B in the (complex) Jordan form B = Bb(BB*, such that MbP = PB, where the (non-singular) 
matrix P has the form 


P = 


Pi -Pi\ 
\P2 pfj- 


This implies that 


{Pi + JIP2) ^Mm{Pi + JIP2) — Bb — Db -\- Nb , 


where Db and N^, respectively denote the diagonal and the nilpotent parts of Bb. 


(18) 


IV-B. C-LINEAR MATRICES 


Let us now consider C-linear transformations. We can associate to any n x n C-linear quaternionic 
matrix a 2n-dimensional complex matrix by the following mapping 


Me = Mb + M;i?, ^ M^ = /(Mb) + I f{M() 


(Ml -M| 
\M2 Mf 


( M'l 


(19) 


where / denotes the isomorphism defined in ( 0 )- Then the following proposition holds: 
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Proposition 1 . 

Let Me be a C-linear quaternionic matrix and its complex counterpart (see Eq. The mapping 

g : Mc<—> Me (20) 

is an isomorphism of the ring of the n-dimensional C-linear matrices into the ring of 2n-dimensional 
complex matrices. 

Indeed, if 

~ ^o,M “P nnd Be — Bq^^ -P 

are two C-linear matrices, their corresponding complex counterparts are given by 

^(^c) = /(^ 0 ,h) + g{Bc) = /(-Bo,h) + ifiBi,a) 

Then, 

-dc Be Ce 

= P1o,h-Bo,h ~ P1i,h-Bi,h + (^0,h-Bi^h + P1-1,h-Bo,h)^2 

and 


g{Ae)g{Be) = /(Ao,h)/(So,h) - /(Ai,h)/(Si,h) -P if{Ai^s)f{Bo,M) + */(^o.h)/(Si,h) 

~ f {-^0 ,bBo^M — Al^gBl^g) -p ifl^AQ gBl g “p Ai^hBq^h) 

= ff(C'c)n 

By using this isomorphism, the right eigenvalue spectrum of C-linear quaternionic matrices can easily 
be determined j^]. The following result 

“A C-linear matrix is diagonalizable if the corresponding complex counterpart is diagonalizable”. 

was proven in Ref. |Q| (where a preliminary discussion of the eigenvalue problem for C-linear quater¬ 
nionic matrix operators was given). It is worth pointing out that the converse of the previous statement 
is, in general, not true. For instance, let us consider the complex matrix 


/zi 0 1 0 \ 

0 Z2 0 1 

0 0 01 0 

\ 0 0 0 02 / 


This matrix admits a corresponding diagonalizable C-linear quaternionic matrix given by 


Ge — De -p Ne — 


/Re( 0 i)-P lm( 0 i)i?i 0 

I 0 Re(02)-P Im(02)i?i 


+ 5 


—j -P kRi 0 
0 -j -P kRi 


( 21 ) 


( 22 ) 


where Ne is nilpotent, diagonal and commutes with De- 

Che normal form of a C-linear quaternionic matrix can easily be calculated. Indeed, given any 
C-linear transformation Me and its corresponding complex counterpart Me [see Eq. (H)] , from the 
known properties of the Jordan form of complex matrices, we can immediately obtain 

S-^MeSe = Je = be + Ne, 


where De is diagonal, Ne is nilpotent and [De, iVc] = 0. Then, the quaternionic C-linear matrices 
Me, Se, De, and Ne are uniquely determined by the isomorphism stated in Proposition 1. 


IV-C. R-LINEAR MATRICES 

In the n dimensional quaternionic vector space P[n,]HI], the R-linear transformations are represented 

by 

3 

Mg_ = ^ H Rfj^ , 

fi—0 


(23) 
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where represent H linear quaternionic matrices and are the right acting operators defined in 
the second section. Any is then characterized by 16n^ real parameters. We can translate M-linear 
nxn matrices into equivalent 4n x 4n real matrices, and vice-versa, by the following translation rules: 


R^ 

Rj 


Rk 


/ 0 - 1 „ 0 0 \ 
1„ 0 0 0 
0 0 0 1 „ 
Vo 0 -1„ 0 / 

/ 0 0 -1„ 0 \ 
0 0 0 -1„ 

1„ 0 0 0 
V 0 1„ 0 0 / 

/ 0 0 0 -1„\ 
0 0 1„ 0 
0 -1„ 0 0 
Vi„ 0 0 0 / 


(24) 


and 


— -^0 ~\~ j ^2 ^ — 


/Mg 

-Ml 

— M 2 

-Mg 

Ml 

Mg 

-Mg 

M 2 

M 2 

Mg 

Mg 

-Ml 

VM3 

-M 2 

Ml 

Mg 


(25) 


where G M[n,H], Mo ,...,3 € M[n,]R], and G M[4n,R]. It is easy to verify that I,J,K 

commute with 


I) — J , 


K, M„ 


= 0 


and 


l2 = J2 = = KJI = -1. 


The following proposition holds. 

Proposition 2 . 

Let Mr be a M-linear matrix and Mr its real counterpart; then the mapping 

h : Mr = Mo,h + Mi^fiRi -P M2^tiRj + M^^^Rk '—> Mr = Mo,h + .^i,hI + LI2,b^ + M3^nK 

is an isomorphism of the ring of the n-dimensional W-linear matrices into the ring of An-dimensional 
real matrices Mr. 

Observe that Let 

Ar = “P -p jA 2 ^wRj T and 7 ?r = .So,h T Ri^wRi B 2 ,uRj T R^^wRk 

be two M-linear quaternionic matrices. Their corresponding real counterparts are given by 

h{^^^ = Ag^e “P IAi,h T JA 2 ^]h[ “P KAg^jj and hl^Bif) = Bq^b -p -p 3B2^m T K7?3^ri 


^rBb - Cr 


= (Ao^niJg^H — Ai^gBi^fi — A2 ^bB2,m Ag^HSg^ii) -p 
(Ag^]H[i?l^e -p Ai^Rji^Q^R -p A2^bB^^b A^^BB2,M)Ri “P 

(Ao^]H[i?2,H “P A2^Mi?0^H “P Ag Al^]H[i?3^]H[)i?j -p 

(Ao_h.B3_h -P Ag^aBo.H + Ai_Hi?2,H — A2,H-Bl,H).Rfe 


Then, 



Stefano De Leo et al: QUATERNIONIC EIGENVALUE PROBLEM 


and 


h{Ag)h{B^) 


(^0, 


- 


- 

^2,h 

-^2,h 

“ A3^h/?3,h) + 

(^0, 


-p 


-p 

^2,h 

-^3,h - 

“ A3^h.B2,h)I + 

(^0, 

h^2,m 

-p 

^2,h-^0,h 

-p 

^3,h 


“ Ai_Hi?3_H)J + 

(^0. 

h^3,m 

-p 

^3,h-So,h 

-p 


^2,h " 

~ A2,h.Si,h)K 

h{C 

ir)^ 








We now discuss the canonical forms of R-linear matrices. Let A* be an R-linear transformation, its 
real counterpart, {Ai + ifii, X 2 + Xs + the complex eigenvalues of A*, and {A 2 S+ 1 , A 4 „} 

the real eigenvalues of A,^. As well known 0,0, there exists a real orthogonal matrix O such that 


where 



("A. 

\ 

P 



Jr — 


X, 




0 

A 2 S +1 




1 

A4n j 




("A. 

\ 

0 


0 

\ 

p 

_ 


A, 

-P 


0 


0 

A 2 S +1 


0 

0 


1 

A4n j 


V 

0 / 


— . 


(26) 


In the previous equation, Xr represents a 2 x 2 real matrix with eigenvalues A,. ± i^r- An appropriate 
choice of O guarantees that 

Let us come back to the R-linear transformation A*. By using the translation rules given in Eqs. 
we can immediately give its canonical form 

Js. — “P Ar . 


In particular, the diagonal elements of Dr corresponding to the quaternionic translation of the real 
blocks 




D 2 , 


D 3 , 


( Ayn \^Ta 0 0 

Mm Ayyi 0 0 

0 0 Xm+1 Mm-t -1 

\ 0 0 Mm+1 Xm+1 

( Am Mm 0 ^ \ 

Mm Am 0 ^1 

0 0 Am+1 0 

V 0 0 0 Am+2 / 

/Am 0 0 0 \ 

0 Am+1 0 0 

0 0 Am+2 0 

V 0 0 0 Am+3/ 



(28) 
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are respectively given by 

-Di,r = ^ [Am(l ~ LiRi) + IIjniLi + Ri) + 

Am+l(l + LiRi) + l-lm+l{Li — Ri)] , 

R2,r — 2 [Am(l LiRi) “P ^m)Li -p Ri) -p 

5 Am+i(l + LiRi — LjRj + LkRk) + 5 Am+2(1 + LiRi + LjRj — L^Rk)] , (29) 

L^3,s. ^ [Ay7i(l LiRi LkRk) ~P A^_|_]^(l LiRi ~P ~P LkRk) “P 

Am+2(1 + LiLli — LjRj + LkRk) + Am+3(1 + LiRi + LjRj — LkRk)] ■ 

As happens for C-linear quaternionic matrices, an R-linear quaternionic matrix is diagonalizable if 
the corresponding real counterpart is diagonalizable. The converse is not necessarily true. 


V. THE EIGENVALUE PROBLEM FOR R-LINEAR MATRICES 

Let us now consider the eigenvalue problem for R-linear quaternionic matrices. Eq. ( |l^ is obviously 
too restrictive. In fact, such an equation sets the real eigenvalue spectrum of R-linear quaternionic 
operators. No information is given about the remaining eigenvalues. In particular, if the real coun¬ 
terpart Mr of the R-linear quaternionic matrix Mr does not have real eigenvalues, Eq. (Hi does not 
admit solution. This is very embarrassing if we consider, for example, R-linear anti-hermitian quater¬ 
nionic operators. Thus, we need to modify Eq. (HD- The discussion regarding the “pseudo-triangular” 
form of the matrices Ar [see Eq. (^^] suggests as R-linear eigenvalue problem the following system 
of coupled equations 

Mg-tp = aij}+ b(p 

Mgip = cip + dip ^ 


where 


a, 6, c, deR, ip = ipo + iipi + jip 2 + kip3 , cp = tpo + itpi + jip 2 + kips , 'i/'o,.... 3 , p5o,...,3 G E[n,R] . 

It can be shown that the real coefficients a, 5, c, d are related to the real and imaginary part of the 
Mr eigenvalues. In fact, by translating the system (B® into its real matrix counterpart, we find 


Mr 0 
0 Mr 


- f bl 




4n 


d4n ^l4n 


(31) 


where 


iP = 


/M\ 


/^o\ 

ipi 


V’l 

i>2 

Vm/ 

, 

^2 

\>f3/ 


G y[4n,] 


The matrix equation (H) admits non trivial solutions if and only if 


det 


ALr 0 , 14 ^ ^l 4 n \ 

ATr d\^ji J 


= 0 . 


(32) 


By rewriting the matrix Mr in terms of the symilarity matrix S'r and of its Jordan form Jr, i.e. 

Mr = Sjj-^ , 

and by using the cyclic property of the determinant, we reduce Eq.(^^ to 


det 


Jr ^l4n 

Jr d^4n J 


1 


= 0 . 


(33) 
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By simple algebraic manipulations |[^, we obtain 


det 


t/ft 

cl47i d^4.n 


— det[(t/R nl47i)((/iK dl47^) 60147^,] 

= det[Jg — (a + d)Jg + {ad — &c)l 4 „] 
= W[z1 - (a + d)zi + {ad - be)] 

i 

= 0 , 


where Zi represent the eigenvalues of the real matrix Mr. The previous equation explicitly shows 
the relation between the real coefficients a, 5, c, d (which appear in the R-linear eigenvalue problem) 
and the eigenvalues z of the real counterpart of the quaternionic matrix M*. In the case of complex 
eigenvalues z, we find 

A = {a + d) — 4(ad — be) = {a — d)^ + be < 0 . (34) 

This condition guarantees that the eigenvalues of the real matrix 

^ _ I al477 bl4n 

y CI 477 dl4n 


appear in conjugate pairs. Consequently, we can find a real similarity transformation T such that 

^ ^ _ I M^4n 

y /4l4n '^l4n 

Finally, without loss of generality, we can consider the following eigenvalue problem for M-linear 
transformations 



MrV' = a '0 - m I 

= a (/? + /r 0 . 


(35) 


VI. FINAL REMARKS 

These final remarks aim to give a concluding discussion on the “coupled” eigenvalue problem and 
a brief summary of mathematical and physical applications motivating our interest in this research. 
In particular, we are interested to bring together two areas: quaternionic differential operators and 
quantum mechanics. 


Coupled eigenvalue equations 


In the previous Section, we have introduced, for R-linear transformations, the eigenvalue problem (35) 
which represents the natural generalization of dil)- In particular, as we observed above, the study 
of system ( |^ instead of Eq. (H) allows to take into account the existence of complex eigenvalues 
and, consequently, complete the eigenvalue spectrum of R-linear quaternionic operators. Actually, the 
eigenvalue problem (|^) also applies to H and C-linear transformations. It can be considered as an 
equivalent formulation of Eqs.(^ and dl)- To show that, let us consider the equation 


Metjj = 'ijjz = V'A + 


(36) 


We limit ourselves to discuss C-linear transformations. Obviously, if a preferred complex direction is 
chosen for the eigenvalues of H-linear quaternionic operators, all the arguments in what follows also 
hold for H-linear transformations. By using the C-linearity, we find 


{Mcip)i = Mc{tpi) = V'jA — 0/4 . 

The pair of eigenvector (0, (^ = 0z), where 0 is solution of Eq.(p6|), satisfies the coupled equations 


Mc0 = A 0 — /r (/? , 

Meyo = /i0 + A(/3 . 


(37) 
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Vice-versa, let Me and be respectively a C-linear transformation and the corresponding eigenvalue 

problem. We denote by <p) a solution of the system If ^ satisfies Eq.(^) too, then, comparing 
Eq.(^^ and the first equation in (|^), one immediately obtains ip = ipi. If, on the contrary, ij) is not 
a solution of Eg. (|3^) , by using the C-linearity, we obtain 

Mc{ipi) = [Mcp)i = pipi + ■ 


Thus, 


Mc('0 -P pi) = ('0 -P pi)z. 


(38) 


Hence, it is possible to associate to any solution {tjj, p) of the system ( |37[) a corresponding eigenvector 
of Me- 

It is worth pointing out that the coupled system ( p5| ) can be obtained by solving the eigenvalue 
problem jl^ ) for complexified quaternionic eigenvectors and complexified real eigenvalues. Infact, by 
imposing that 

tjj ^+1p e ]Hl(i, i, j, fc) (g) C(i,/) , 

and 

r ^ Z = X + Ip G C(1,J) , 
from the complexified eigenvalue problem (H), 


we immediately get the coupled system (35). 


Applications 

Many physical problems dealing with differential operators are greatly simplified by using the matrix 
formalism and solving the corresponding eigenvalue problem. 

Let us first consider a very simple case. That is the H-linear second order homogeneous ordinary 
differential equation 


'fix) — a if{x) — /3 'fix) = 0, q;,/3g]HI,xGK.. (39) 

In looking for quaternionic exponential solution tfix) = exp[qx] and observing that the derivative of 
exp[ga;] with respect to the real variable x is q exp[( 7 a:], we reduce the previous problem to find the 
solutions of the following quadratic equation 

= aq + (3 . (40) 

This equation can be rewritten in matrix form as follows 



As seen in this paper, the H-linear quaternionic matrix satisfies a right (complex) eigenvalue 
equation 

= z eC , v,w eU , wiM = 1 . (42) 

Due to the particular form of M^, the components of the Mn-eigenvectors satisfy the following condi¬ 
tion 

V = cu z . (43) 

Multiplying (from the right) Eq. (^^ by re and using the constraint (^), we obtain 




/ wzw 

V 1 


I WZ^IV 
WZ'W 


(44) 
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Comparing Eq. (El with Eq. (Q), we immediately get 

q = w zw . 


( 45 ) 


The problem of finding exponential solutions for H-linear differential with constant coefficients and, 
consequently, zeros of H-linear polynomial equations ||^ , is thus equivalent to solve the right (com¬ 
plex) eigenvalue problem for the associated matrix. Obviously, the previous considerations also hold 
for the n-dimensional case. 

The solutions of X-linear quaternionic differential equations with constant coefficients 

(x) - (x) - (x) - ... - Ao,,V'(^) = 0 , X = R, C, H . (46) 


can be given in terms of the eigenvalues and eigenvectors of the matrix 


/A 


n—l,x A^_2 ,x 

1 0 


Ao,x ^ 
0 


V 0 


0 . . 1 


0 


Interesting C-linear differential equations appear in quaternionic quantum mechanics . Eor example, 
by studying quaternionic tunneling effect as candidate to possible phenomenological deviations from 
the standard (complex) theory, we have to solve the following H-linear Schrodinger equation 


9t4'(x,f) 






4'(x,t) 


(47) 


where j- W represents the new quaternionic perturbation. The quaternionic stationary state wave 
function 

4>(x, t) = ^/>(x) exp [—I Et] (48) 

is solution of Eq. ( ^^ on the condition that ipix) be solution of the following time-independent C-linear 
(ordinary) differential equation 

i ^ '^(x) — iV V'(x) + jW ^(x) -P ipix) i E = 0 . (49) 

Observe that the choice of the imaginary unit i in the Laplacian operator Eq. (p^ , and in the 

time exponential, Eq. (^^ , is fundamental to recover the standard results in the complex limit. In this 
formalism, quaternionic potentials are treated as perturbation effects on standard quantum mechanics. 
We also point out that the right position of the time exponential is fundamental to perform the 
separation of variables. 

The solution of Eq.(^^ can be given in terms of the eigenvalues and eigenvectors of the C-linear 
matrix 

fOAoA 

Vi 0 y ’ 

where 

^o,c = ^ {V + LkW + L,R,E) . 

A detailed phenomenological discussion of the quaternionic tunneling effect is found in the paper of 
ref. i. 


Outlooks 

As seen in this paper, the choice of right (complex) eigenvalues for H and C linear operators play 
a fundamental role in discussing canonical forms, in finding solutions of polynomial and differential 
equations. It was shown that the right (complex) eigenvalue problem is equivalent to a “coupled” 
system and this was extremely important to study the eigenvalue problem for M-linear quaternionic 
matrices, where a pair of real eigenvalues must be introduced. This work was intended as an attempt 
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at motivating the study of R and C-linear quaternionic operators in view of possible applications in 
quantum mechanics and gauge theory. It would be desirable to give a complete theory of X-linear 
quaternionic matrices and differential operators. More realistically, this paper touches only a few 
aspects of the theory and shows how the choice of the right eigenvalue equation seems to be the 
best to investigate quaternionic formulations of physical theory. It was not our purpose to study 
here differential operators. The results in this field are far from being conclusive and some questions 
represent at present intriguing challenges: variations of parameters; order reduction; not invertible 
higher derivative K and C-linear constant coefficients; variable coefficients; integral transforms. Finally, 
it would be desirable to extend the discussion on the eigenvalue problem by matrix translation to the 
non-associative case 
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